Modeling bed erosion in free surface flows by the particle finite element method

We present a general formulation for modeling bed erosion in free surface flows using the particle finite element method (PFEM). The key feature of the PFEM is the use of an updated Lagrangian description to model the motion of nodes (particles) in domains containing fluid and solid subdomains. Nodes are viewed as material points (called particles) which can freely move and even separate from the fluid and solid subdomains representing, for instance, the effect of water drops or soil/rock particles. A mesh connects the nodes defining the discretized domain in the fluid and solid regions where the governing equations, expressed in an integral form, are solved as in the standard FEM. The necessary stabilization for dealing with the incompressibility of the fluid is introduced via the finite calculus (FIC) method. An incremental iterative scheme for the solution of the nonlinear transient coupled fluid-structure problem is described. The erosion mechanism is modeled by releasing the material adjacent to the bed surface according to the frictional work generated by the fluid shear stresses. The released bed material is subsequently transported by the fluid flow. Examples of application of the PFEM to solve a number of bed erosion problems involving large motions of the free surface and splashing of waves are presented

Modeling bed erosion in free surface flows by the particle finite element method

We present a general formulation for modeling bed erosion in free surface flows using the particle finite element method (PFEM). The key feature of the PFEM is the use of an updated Lagrangian description to model the motion of nodes (particles) in domains containing fluid and solid subdomains. Nodes are viewed as material points (called particles) which can freely move and even separate from the fluid and solid subdomains representing, for instance, the effect of water drops or soil/rock particles. A mesh connects the nodes defining the discretized domain in the fluid and solid regions where the governing equations, expressed in an integral form, are solved as in the standard FEM. The necessary stabilization for dealing with the incompressibility of the fluid is introduced via the finite calculus (FIC) method. An incremental iterative scheme for the solution of the nonlinear transient coupled fluid-structure problem is described. The erosion mechanism is modeled by releasing the material adjacent to the bed surface according to the frictional work generated by the fluid shear stresses. The released bed material is subsequently transported by the fluid flow. Examples of application of the PFEM to solve a number of bed erosion problems involving large motions of the free surface and splashing of waves are presented

PFEM application in fluid structure interaction problems

In the current paper the Particle Finite Element Method (PFEM), an inno-vative numerical method for solving a wide spectrum of problems involving the interaction of fluid and structures, is briefly presented. Many examples of the use of the PFEM with GiD support are shown. GiD framework provides a useful pre and post processor for the specific features of the method. Its advantages and shortcomings are pointed out in the present wor

Accurate modelling of the elastic behavior a continuum with the Discrete Element Method

The Discrete Element Method (DEM) has been used for modeling continua, like concrete or rocks. However, it requires a big calibration effort, even to capture just the linear elastic behavior of a continuum modelled via the classical force-displacement relationships at the contact interfaces between particles. In this work we propose a new way for computing the contact forces between discrete particles. The newly proposed forces take into account the surroundings of the contact, not just the contact itself. This brings in the missing terms that provide an accurate approximation to an elastic continuum, and avoids calibration of the DEM parameters for the purely linear elastic range

A FEM-DEM technique for studying the motion of particles in non-Newtonian fluids. Application to the transport of drill cuttings in wellbores

We present a procedure for coupling the finite element method (FEM) and the discrete element method (DEM) for analysis of the motion of particles in non-Newtonian fluids. Particles are assumed to be spherical and immersed in the fluid mesh. A new method for computing the drag force on the particles in a non-Newtonian fluid is presented. A drag force correction for non-spherical particles is proposed. The FEM-DEM coupling procedure is explained for Eulerian and Lagrangian flows and the basic expressions of the discretized solution algorithm are given. The usefulness of the FEM-DEM technique is demonstrated in its application to the transport of drill cuttings in wellbores

Advances in the particle finite element method for fluid-structure interaction problems

We present a general formulation for analysis of fluid-structure interaction problems using the particle finite element method (PFEM). The key feature of the PFEM is the use of a Lagrangian description to model the motion of nodes (particles) in both the fluid and the structure domains. Nodes are thus viewed as particles which can freely move and even separate from the main analysis domain representing, for instance, the effect of water drops. A mesh connects the nodes defining the discretized domain where the governing equations, expressed in an integral from, are solved as in the standard FEM. The necessary stabilization for dealing with the incompressibility of the fluid is introduced via the finite calculus (FIC) method. A fractional step scheme for the transient coupled fluid-structure solution is described. Examples of application of the PFEM to solve a number of fluid-structure interaction problems involving large motions of the free surface and splashing of waves are presented

Possibilities of the particle finite element method for fluid–soil–structure interaction problems

We present some developments in the particle finite element method (PFEM) for analysis of complex coupled problems in mechanics involving fluid–soil–structure interaction (FSSI). The PFEM uses an updated Lagrangian description to model the motion of nodes (particles) in both the fluid and the solid domains (the later including soil/rock and structures). A mesh connects the particles (nodes) defining the discretized domain where the governing equations for each of the constituent materials are solved as in the standard FEM. The stabilization for dealing with an incompressibility continuum is introduced via the finite calculus method. An incremental iterative scheme for the solution of the non linear transient coupled FSSI problem is described. The procedure to model frictional contact conditions and material erosion at fluid–solid and solid–solid interfaces is described. We present several examples of application of the PFEM to solve FSSI problems such as the motion of rocks by water streams, the erosion of a river bed adjacent to a bridge foundation, the stability of breakwaters and constructions sea waves and the study of landslides

Numerical modelling with discrete elements of rockfall protection systems

Some important infrastructures like roads, railway tracks or dams were constructed in places threatened by natural hazards. With the purpose of preserving these infrastructures from landslides and rock-falls, different containment systems are installed, and one of the most popular are the flexible metallic fences. The development of full-scale laboratory tests to evaluate the behaviour of flexible metallic fences is unfunctional, accounting to the huge magnitude of the event. On the other hand, small-scale testing may lead to inaccurate results, due to the distortion in the contours (e.g. anchors of the metallic fences). These problems in laboratory testing have led to the popularization of the use of numerical methods. In this study, the bonded Discrete Element Method (DEM) is used for the analysis of the behaviour of flexible metallic fences for rockfall protection. The bonded DEM is a modification of the classical DEM which assumes that bonds exist between particles, resisting their separation. In this case, the net cables are represented using rigid spheres joined by bond elements that are deformed according to an elasto-plastic law. Calculations were carried out using the DEMPack program, a specific software developed in CIMNE for modelling with the bonded DEM. This software allows considering the inter-action between discrete and finite elements, which can be useful to represent the boundaries of the domain, such as the surface of the slope. The code is firstly validated reproducing benchmark tests available in the literature. Finally, full-scale tests are computed in order to evaluate the energy dissipation capacity of the fence during a rockfall event

Lagrangian analysis of multiscale particulate flows with the particle finite element method

We present a Lagrangian numerical technique for the analysis of flows incorporating physical particles of different sizes. The numerical approach is based on the Particle Finite Element Method (PFEM) which blends concepts from particle-based techniques and the FEM. The basis of the Lagrangian formulation for particulate flows and the procedure for modelling the motion of small and large particles that are submerged in the fluid are described in detail. The numerical technique for analysis of this type of multiscale particulate flows using a stabilized mixed velocity-pressure formulation and the PFEM is also presented. Examples of application of the PFEM to several particulate flows problems are given